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I. INTRODUCTION 



The R- matrix theory of reactions has proven over the course of time to be very useful in nuclear and atomic physics, 
both for the fitting of experimental data and as a tool for theoretical calculations. In these notes I will explore the 
' relationship between formal 12-matrix parameters and "physical" resonance parameters. 
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CN 1 II. HOW DOES ONE DEFINE A RESONANCE ENERGY? 

X> \ 

In practice one is likely to encounter many definitions of "the resonance energy" - e.g. the energy where there is a 
peak of an excitation function, the energy when the phase shift is equal to n/2 (for elastic scattering), the complex 
OO poles of the scattering matrix,. . .Many approaches involve the Breit-Wigner resonance formula: 

cm : 

T — | ; cc ' 1 ' ~ B c (E R - Ef + (r°/2) 2 { ' 

where c(c') label the incoming (outgoing) channel, k 2 c is the incoming wavenumber, u> = (2J+l)/[(2Ji + l)(2J2 + l)] is 
the statistical factor, E is the (variable) reaction energy, Er is the resonance energy, T° are the observed partial widths 
of the resonance, and the total observed width r° is given by sum over all channels of the partial widths: T° = ^ c T°. 
It should be noted that the Breit-Wigner formula is a very general quantum-mechanical result with applications 
, ranging from atomic physics to particle physics (e.g. the 92-GeV resonance in e + — e~ scattering known as the Zq 
boson). In some cases the partial widths may be energy dependent. From a time-dependent viewpoint the resonant 
state decays with time-dependence exp(— t/r), where r is the mean lifetime and T°t = h. I will describe below one 
■ approach for defining resonance energies and partial widths in the i?-matrix formalism which is particularly convenient. 
The differences between various prescriptions for defining resonance energies and partial widths are generally most 
significant for broad resonances. Let me conclude this section with a caution: since different methods are in use for 
defining resonance parameters, it is very important to (1) understand what conventions are used by others if you use 
their results and (2) supply enough information so that others may properly use your results. 
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III. REVIEW OF .R-MATRIX THEORY 

I will begin by reviewing some of the notation and results of i?-matrix theory as described by Lane and Thomas 
(LT) [1] which will be needed later. The R matrix is a function of the energy E and is defined by 

ij cV = y^ (2) 

^ Ex -E' y ' 

where E\ are the level energies, 7> c are the reduced width amplitudes, A is the level label, and c is the channel label. 
I will assume that the numbers of levels and channels are finite and given by N\ and N c , respectively. One must also 
specify the constant(s) B c , which determine the boundary conditions satisfied by the underlying eigenfunctions. I will 
also consider here a set of levels which all have the same value (the presence of additional levels with other J™ 
values would not impact any results discussed here). 

When working in the level-matrix framework, it is convenient to define column vectors 7 C in level space from the 
7a c : 



T« = [ ~' 12c I (3) 
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as well as a diagonal matrix e containing the Ey. 

[e] AjU = ExS Xfi . 

The N\ x N\ level matrix A can now be defined by its inverse: 

C 

A" 1 = e-El-J2lc7l(S c + iP c -B c ), 



(4) 

(5) 
(6) 



where 1 is the unit matrix and P c and S c are the energy-dependent penetration and shift functions. Note also the 
notation X T = transpose(X) for any matrix X. The collision matrix U which determines the cross sections can be 
calculated from the A matrix: 



[U] 



C'C "C L c 



8 dc + 2i(P c ,P c ) 1 / 2 1 J,A lc 



(7) 



where the fi c are defined in terms of Coulomb wavefunctions (see LT, Eq. III. 4. 5). 



IV. THE SINGLE-LEVEL fl-MATRIX 

In order to interpret the meaning of the _R-matrix parameters E\ and 7 Ac it is very helpful to consider the case of 
a single level, i.e. N\ = 1. In this case the A matrix is simply a number and the cross section is given by 



<7 CC >(E) = 



TxrT 



Ac^ Ac' 



kl (s A + A A -i?)2 + (r A /2)2' 



where 



r Ac - 2jI c p ( 

A A = 



r A 



-5>AcPc(£)-Bc], and 



(8) 



(9) 
(10) 

(11) 



The partial and total widths r Ac and T A are known as formal widths. This formula for the cross section is very 
similar to the Breit-Wigncr formula defined above by Eq. (1) (in fact the single-level U-matrix formula is sometimes 
called the Breit-Wigncr formula). The relationship between Eqs. (1) and (8) is simplest if the B c are chosen such 
that B c = S C (E\). The energy-dependent quantity A A can then be expanded in a Taylor series around E — Ey. 



(E-Ex 



Ac 
dE 



+ ... 



dS r 



If we replace A A with the first non-zero term in the expansion, we have 



(12) 
(13) 



(E x + A x - Ef » (E x - Ef 



2 i dS c \ 

^Ie) 



We can now connect the i?-matrix parameters with the Breit-Wigner resonance parameters. If we identify 

2 1 l c P c (Ex) 



E R = E x and 

r Ac (£ A ) 



r? 



1 + Ec7L(^) £a 1 + E C 7L(^) 



(14) 

(15) 
(16) 



3 



then the one-level _R-matrix formula [Eq. (8)] and the Breit-Wigner formula [Eq. (1)] are equivalent. The presence of 
the penetration factor in Eq. (9) gives rise to an energy dependence of the I\ c hence the need for defining F\ C (E\) 
and P C (E\) in Eq. (16). Eqs. (15) and (16) can be considered as the definition of the resonance energy and partial 
widths for the case of a single- level R- matrix (i.e. the physical parameters). These relationships are also discussed in 
Sec. XII. 3(a) of LT. The resonance energy defined in this manner may differ from the peak of cross section excitation 
function due to the energy dependences of the penetration factors and k 2 c . Note also that the difference between 
formal and observed partial widths is often on the order of 30% - a large factor compared to typical experimental 
uncertainties. 



V. GENERAL DEFINITION OF PHYSICAL RESONANCE PARAMETERS 



When there is more than one level, a one-level approximation to the R matrix is very good when E is near 
the R- matrix pole energy E\ of a particular level [this approximation becomes exact when E — > Ex; see Eq. (2)]. 
Eqs. (15) and (16) from the previous section thus define the physical parameters for that particular level. However, 
this approach requires defining the boundary condition constants such that B c = S C (E\). The B c can only be set 
once (they are after all "constants" ) - therefore we only have a simple relationship between i?-matrix parameters and 
physical parameters for one of the levels. 

One approach to finding the physical parameters for the remaining levels is to change the B c to values appropriate 
for another physical resonance energy. The _R-matrix parameters E\ and jxc must then be re-determined - either 
by refitting data or by the transformation procedure discussed in the next section. This procedure must be then 
be repeated one level at a time until physical parameters have been found for all of the levels. Note also that this 
procedure is necessarily iterative - we want to set B c — S C {E\) where E\ is a physical resonance energy, but the 
physical resonance energy is one of the quantities we are seeking to determine. Fortunately an iterative approach 
starting with the H-matrix pole energy converges quickly. 



VI. HOW TO CHANGE THE B c 



An interesting feature of i?-matrix theory is that the collision matrix (i.e. the U matrix) is invariant under changes 
in the B c , provided that the E\ and 7a c are suitably adjusted. This result remains true even for the case of finite 
N\ [2]. The transformation is most easily described using matrix equations in level space. Let us consider the 
transformation B c — > B' c , E\ — > E' x , and -jxc — ► j' Xc . One first constructs the real and symmetric matrix C defined by 

C = e-Y,JclI(B , c -B c ), (17) 

c 

which is diagonalized by the orthogonal matrix K such that D = KCK J ', with — D\5\^. Note that any real 
and symmetric matrix can be diagonalized by an orthogonal matrix; an orthogonal matrix K satisfies K T = K" 1 
(in other words K T K = KK T = 1). The necessary transformation of the i?-matrix parameters is then given by [2] 

E'x = D X (18) 

and 

7c - K lc . (19) 

By using Eq. (6) one can show that these transformations leave ^A-jc and hence U [see Eq. (7)] invariant. 

The invariance of U for changes in B c has several important consequences. I will highlight two of them: (1) 
mathematically-equivalent cross section fits should be attainable for any B c choices and (2) we are free to choose B c 
however we desire. We can exploit this freedom by adopting B c values which are convenient [e.g. B c = S C (E\)]. This 
situation is somewhat analogous to that encountered in electromagnctism: transformations of the scalar and vector 
potentials (known as gauge transformations) can be carried out which leave the physical electric and magnetic fields 
invariant. 

It is instructive to look at the transformation for the case of a single level. The matrix C consists of only a single 
dimension, i.e. a number. The resulting transformation is 

E x = E.-^UB'c-Bc) (20) 

C 

IXc = 7Ac (21) 
It is straightforward to verify by substitution into the single-level formula [Eq. (8)] that the cross section is invariant. 
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VII. ANOTHER APPROACH TO PHYSICAL PARAMETERS 



In this section I will describe another mathematically-equivalent method for extracting physical parameters from 
H-matrix parameters. I begin by defining the real and symmetric matrix £: 

S = e-Y / lclI[Sc(E)-B c }, (22) 

c 

and consider the eigenvalue equation 

Sa, = Eitti (23) 

where Ei is the eigenvalue and ai is the corresponding eigenvector. Note that £ is implicitly dependent upon 
Ei through S c , so the eigenvalue problem is nonlinear. We will assume for convenience that the eigenvectors are 
normalized so that afdi = 1. 

The eigenvalues Ei are in fact just the physical resonance energies. A new set of reduced width parameters is 
defined via 

il c = af lc . (24) 

The 7i C are equal to the i?-matrix reduced widths when transformed to B c = S c {Ei). The ji C are sometimes called 
the "on-resonance" reduced-width amplitudes. The observed (physical) partial widths are thus given by 

2P c (E l )jf c 

1 ^ Z^c he \ dE ) Ei 

One advantage of this formulation is that the solution of the eigenvalue problem gives immediately all of the physical 
resonance energies and partial widths. 

However one defines the physical resonance parameters, it is desirable that they be independent of the boundary 
condition constants and channel radii of i?-matrix theory. The above approach is completely independent of the B c 
and approximately independent of the channel radii. Additional details and the justification for this approach are 
given by Brune [3] . The "on-resonance" reduced- width amplitudes defined here are also disused by Barker [4] and by 
Angulo and Descouvemont [5]. 

I would like to mention in passing another similar eigenvalue equation. If we had instead defined 

£ = e - ^ 7c7 J[S c (£) + \P C (E) - B c ], (26) 

c 

we would have the eigenvalue equation introduced by Hale, Brown, and Jarmie [6]. The eigenvalues of this equation 
are complex and correspond to complex poles of the A and U matrices (the latter is also known as the scattering 
matrix). As discussed in Refs. [6, 7] this approach to defining resonance energies and widths may be theoretically 
preferable. However, there are significant complexities associated with calculating the Coulomb wavefunctions for 
complex momenta. We will not discuss this approach further here. 



VIII. CONVERTING PHYSICAL PARAMETERS TO .R-MATRIX PARAMETERS 

Next we will investigate how to translate physical parameters into i?-matrix parameters [3] . It is often desirable to 
incorporate information from a level compilation or shell model calculation into an i?-matrix parameterization. Values 
for r° c must first be converted into reduced widths ji C using Eq. (25). This procedure requires the straightforward 
solution of a system of linear equations for jf c ; the only complication is that the signs of reduced- width amplitudes 
cannot be determined. Relative signs can be determined via an H-matrix calculation (which automatically includes 
interference effects between levels) and a comparison to experimental data. 

With the physical resonance energies Ei and "on-resonance" reduced- width amplitudes 7i C in hand we can determine 
i?-matrix parameters. We proceed by defining the matrices M and N: 

Mii = \ 1 y z. z. s ic -s jc \~J (27) 
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and 

r^i I Ei + Ec^USic-Bc) i = j 

[N]ij = \E c i ic 7 jc [^E^-B c ) i*j> (28) 

where S ic denotes S c (Ei). Note that the construction of AT requires the adoption of specific B c values. The eigenvalue 
equation 

Nb x = E\Mb\. (29) 

now holds the key for transforming from the Ei-j ic representation to the standard i?-matrix parameters E\ and 
7a c . The eigenvalues are the R- matrix pole energies E\ and the eigenvectors 6a yield the il-matrix reduced- width 
amplitudes: 

7Ac - bl%. (30) 

The matrices M and TV are real, symmetric, and energy-independent. The solution to this type of eigenvalue problem 
is discussed in Sec. 8.7.2 of Ref. [8]; I have utilized the LAPACK [9] routine dsygv. 

The eigenvectors b\ can be arranged into a square matrix b which simultaneously diagonalizes M and N: 

b T Mb = 1 (31) 

and 

b T Nb = e. (32) 

We can also re-write Eq. (30) as 

7c = b T %. (33) 



IX. WORKING DIRECTLY WITH PHYSICAL PARAMETERS 



It turns to be possible to carry out R- matrix calculations directly from the physical parameters [3] . This result can 
be achieved by defining a new (alternative) level matrix A from the physical resonance energies E t and "on-resonance" 
reduced- width amplitudes -ji C : 



[A-% = {E i -E)5 ij -Y,liclic{Sc + iPc) 



(34) 



E 



llc S ic * = 3 

lie Ijc Ei-Ej ' J 



This matrix is then used with the 7 C to calculate the the collision matrix U: 



5 c , c + 2i(P c ,P c ) 1 / 2 jJ,Aj c 



(35) 



This method of calculating U is mathematically equivalent to first transforming the physical Ei-"fi C parameters to 
J£-matrix parameters E\-"/\ c and then calculating U using standard _R-matrix equations. Note that in this approach 
the boundary-condition constants B c never need to be defined at all. One application of this approach is the possibility 
of using the physical parameters directly as fit parameters. 



X. AN EXAMPLE 



At this point it is instructive to introduce a concrete example. The top section of Table I shows a simple well- 
documented set of standard il-matrix parameters taken from Azuma et al. [10]. This analysis concerns the low-lying 
,V = 1~ states of 16 reacting through channels 12 C + a (c = a, l a = 1) and 16 + 7 (c = 7, El multipolarity) . The 
zero of the energy scale is at the 12 C + a threshold. Note that the gamma-ray channel is treated perturbatively and 
does enter in the calculation of the A matrix. 
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TABLE I: Standard ii-matrix parameters from Table III of Ref. [10] which describe J n = 1 reactions in the le O system, 
and the physical parameters derived from them. The channel radius is 6.5 fm and the boundary condition constant is chosen 
according to B c = S C (E 1 ). The /3-decay feeding amplitudes Bx are equivalent to the quantities A\i~/^ Na 1,72 of Ref. [10]. 



A 


1 


2 


3 


E x (MeV) 


-0.0451 


2.845 


11.71 


7Aa (MeV 1 / 2 ) 
7a 7 (MeV -1 ) 


0.0793 
8.76 x 10~ 6 


0.330 
—2.44 x 10 


1.017 

—2.82 x 10~ 6 


Bx 


1.194 


0.558 


-0.629 


E x (MeV) 


-0.0451 


2.400 


8.00 


7a« (MeV 1/2 ) 
7 a 7 (MeV" 1 ) 


0.0793 


0.471 


0.912 


8.76 x 10~ 6 


-3.20 x 10~ 6 


-2.50 x 10~ 6 


Bx 


1.194 


0.408 


-0.781 


T° (MeV) 
T° (meV) 


1.0050 
55.1 


1.2034 

0.359 

8.64 


1.0220 
9.041 
54.2 



TABLE II: Elements [b]ij of the transformation matrix corresponding to the parameters of Table I. 



\ 3 
i \ 


1 


2 


3 


1 


1.000 


0.0373 


0.0446 


2 


0.000 


0.9781 


0.2281 


3 


0.000 


-0.1466 


0.9933 



The bottom section of Table I shows physical resonance energies Ei and the "on-resonance" reduced-width ampli- 
tudes 7a c and /3-decay feeding factors Bi. In addition the physical partial widths T° c are given along with values of 
[1 + 7, 2 a ^(Bi)] needed to related formal and physical widths. 

This example is also useful for demonstrating some of the other concepts discussed here. For example one can 
convert from the Ei-^ ic representation back to the standard ii-matrix E\-"f\ c representation. In Table II we show 
the elements of the transformation matrix b [see Eqs. (31-33)] needed to carry the transformation back to standard 
il-matrix parameters. One can also verify that the various approaches for calculating the collision matrix U [e.g. 
Eqs. (7) and (35)] yield the same result. 



XI. SUMMARY 



Several aspects of the relationship between il-matrix parameters and physical resonance parameters have been 
discussed. First it is necessary to establish how physical resonance energies and partial widths will be defined. 
By comparing the one-level _R-matrix formula for the cross section to the Breit-Wigner formula it is possible to 
define a relationship between _R-matrix parameters and physical resonance parameters. For the multi- level case the 
extraction of physical parameters is more complicated but can be achieved by the refitting data with different B c values, 
carrying out B c transformations, or by solving an eigenvalue equation. We next investigated the reverse procedure, i.e. 
transforming physical resonance parameters into il-matrix parameters, which is useful for incorporating information 
from other sources into an il-matrix fit. A method for calculating observables directly from the "on-resonance" 
il-matrix parameters was then discussed. Finally we presented an example case which is useful for demonstrating the 
ideas which have been presented here. 
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